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We define a local Riemannian metric tensor in the manifold of Gaussian channels and the distance 
that it induces. We adopt an information-geometric approach and define a metric derived from the 
Bures-Fisher metric for quantum states. The resulting metric inherits several desirable properties 
from the Bures-Fisher metric and is operationally motivated from distinguishability considerations: 
It serves as an upper bound to the attainable quantum Fisher information for the channel parameters 
using Gaussian states, under generic constraints on the physically available resources. Our approach 
naturally includes the use of entangled Gaussian probe states. We prove that the metric enjoys 
some desirable properties like stability and covariance. As a byproduct, we also obtain some general 
results in Gaussian channel estimation that are the continuous- variable analogs of previously known 
results in finite dimensions. We prove that optimal probe states are always pure and bounded in 
the number of ancillary modes, even in the presence of constraints on the reduced state input in 
the channel. This has experimental and computational implications: It limits the complexity of 
optimal experimental setups for channel estimation and reduces the computational requirements 
for the evaluation of the metric: Indeed, we construct a converging algorithm for its computation. 
We provide explicit formulae for computing the multiparametric quantum Fisher information for 
dissipative channels probed with arbitrary Gaussian states, and provide the optimal observables for 
the estimation of the channel parameters (e.g. bath couplings, squeezing, and temperature). 

PACS numbers: 03.65.Ta, 03.67.Hk, 42.50.Dv 



I. INTRODUCTION 

The theory of quantum channels provides a broad con- 
ceptual and mathematical framework to describe physical 
transformations on quantum states. Progress in quantum 
information technology is bringing long-standing ques- 
tions related to quantum channels to the front line of re- 
search. Topics such as dissipation-assisted quantum com- 
putation fi[ , quantum teleportation 0] , quantum memo- 
ries [U and quantum state engineering [j, Q all have in 
common that they deal with quantum channels in one 
way or another. One major question that has recently 
received significant amount of attention is the definition 
of a distance among quantum channels The main 

motivation for such notion is the identification of a gold 
standard Q against which all quantum processes could 
be compared, which would unify and systematize the way 
in which errors are treated and quantified. In @ a rea- 
sonable, physically motivated set of requirements for such 
distance was introduced and some particular cases fulfill- 
ing most of them were identified. However, despite recent 
progress in the field Q, many questions remain open. 

A similarly motivated program, as developed between 
the late sixties and the early nineties, regarded the prob- 
lem of denning distances among quantum states. The 
resulting theory is today well encompassed within the 
framework of Information Geometry (T^-tHJ. In par- 
ticular, this approach addresses questions regarding the 
distinguishability and the estimation of sets of quantum 



states. These questions led to the notions of quantum 
fidelity [Hj], Bures distance and quantum Fisher infor- 
mation (QFI) HM^. As a result, the manifold of quan- 
tum states is endowed with a local Riemannian metric 
tensor which is physically motivated and serves as the 
gold standard for comparing quantum states. Addition- 
ally, the metric allows one to define the Bures distribu- 
tion fill l20l l2ll ]. the quantum analog to the Jeffreys' 
prior [221 ] . which, among other applications, provides an 
operationally motivated prior distribution for Bayesian 
tomography and estimation techniques [23^ . 

There exists in the literature a number of proposals 
for defining a distance for quantum channels. Most no- 
tably, the Jamiolkowski process distance Q, based on 
the Choi- Jamiolkowski isomorphism [24|, and the com- 
pletely bounded trace norm [J l25j establish a distance 
between two channels Si and £2 by considering the dis- 
tance d(pi,p2) between the respective images of an ar- 
bitrary probe state p under the action of the channels 
Pi = (Si ®T)p. Such distance is then maximized over 
all possible probe states p in order to obtain a funda- 
mental measure of distinguishability between S\ and £2 ■ 
The completely bounded trace norm is arguably the most 
appropriate choice both for physically motivated reasons 
and practical considerations: Indeed, it has been recently 
shown that it can be computed in polynomial time in the 
dimension of the system upon which the channels act 0] . 
Despite this and other significant contributions, existing 
distance measures are inadequate to address a number of 



2 



relevant situations. In particular, defining a distance on 
a manifold instead of a metric tensor hinders the task of 
defining natural prior distributions over the set of chan- 
nels. Moreover, most of the existing distances are hard, 
if not impossible, to compute in infinite-dimensional sys- 
tems such as continuous variables [III [2?J . Additionally, 
in a number of situations, arbitrarily good distinguisha- 
bility between any two infinite-dimensional channels can 
be achieved if one allows for a large enough amount of 
resources [281 ] . thus rendering the direct approach of op- 
timization [9( useless, unless some restrictions on the re- 
sources (regularization) are enforced. 

In the present work we address the problem of defin- 
ing a metric tensor on the set of n-mode Gaussian chan- 
nels £ [ll], and we study the ensuing physical conse- 
quences. The set of Gaussian channels can be regarded 
as a manifold once a proper paramctrization X : O — >• £ 
is established, where O C M. d is an open subset in a 
e?-dimcnsional real vector space. The set £ = {S(X)} 
can be equipped with a Riemannian metric following the 
spirit of the statistical distance [33[ and the Bures dis- 
tance. Considering any two infinitesimally close quantum 
channels <Si = S(X) and S2 = S(X + dX), a metric ten- 
sor 3 provides the infinitesimal distance between them 
as 

#{Sx,S 2 )=Z lu ,{X)dX> l dX v (1) 

where 3{X) > 0. We use Einstein's summation conven- 
tion, i.e., an index appearing once as a sub- and once as 
a super-index is automatically contracted. 

As pointed out in Q some physically motivated re- 
quirements should be imposed. Namely, the metric 
should meet the following criteria: 

• Stability: The metric should be invariant under the 
addition of ancillary modes, namely d 2 {S\,S2) = 
d 2 (Si ®T, S2 where X corresponds to the iden- 
tity channel on an arbitrary number of ancillary 
modes. 

• Measurability: This amounts to say that, once the 
metric is defined, the channel parameters can be 
determined by experimental means. 

• Computability: The metric Z{X) should be com- 
putable. This requirement is obviously unclear 
for the existing proposals when addressing infinite- 
dimensional channels. 

• Physical meaning: The metric should have a clear 
operational meaning. This is certainly the case if 
the metric is derived from distinguishability con- 
siderations. 

Other requirements (symmetry, non-degeneracy, and the 
triangle inequality) are guaranteed by any distance stem- 
ming from a metric tensor. The chaining condition im- 
posed on Q is not immediate to translate into the metric 
approach, and we will not address it in the present work. 
Finally, a purely formal requirement is in order for any 
well-defined metric tensor 



• Covariance: The metric tensor has to transform 
covariantly under a reparametrization. Namely, if 
we perform the change X —> X'(X), with dX'^ = 
A^ u (X)dX v the metric tensor has to transform as 
VPO = 3' X<T (X')A\(X)A° U (X) so that the dis- 
tance d(S±,S2) is invariant under reparametriza- 
tion. 

As mentioned earlier, some physical constraint on the 
probe states needs to be imposed in order to guarantee 
that the metric is well defined. We assume that this 
constraint is given in the form of a real- valued function 
<A(p) < 0* , where the specific form of <fi is not of particular 
relevance. This function may have a practical motivation, 
and thus be chosen according to technical considerations, 
or may be used to regularize the divergencies that appear 
when resources are unlimited; we will refer to it as a 
resource budget. From now on we tacitly assume that all 
quantities we define depend implicitly on the given choice 
of 0. 

The paper is organized as follows. SectionlTIlintroduces 
the Bures-Fisher Metric, a concept central to our work, 
together with some of its properties. Section IIIII defines 
our channel metric and derives some of its main prop- 
erties, namely covariance and stability, as well as show- 
ing that it can be computed with arbitrary precision, 
provided that one can compute Quantum Fisher infor- 
mation (QFI) matrices for all Gaussian probe states. In 
doing so, we generalize some results known previously for 
finite-dimensional channels into the constrained Gaussian 
setup. In Section IIVI we derive explicit formulae for ob- 
taining the QFI matrices for dissipative channels. Section 
fVl provides some remarks and stresses the main ques- 
tions left open, as well as the near-future applications of 
our results to the estimation of relevant channel proper- 
ties such as bath couplings, temperature, and squeezing. 
The mathematical details are reported in five technical 
appendices. 



II. THE BURES-FISHER METRIC 

The main motivation for defining a metric in the 
manifold of Gaussian channels £ is to obtain a 
parametrization-independent measure of distinguishabil- 
ity on £. By expressing distances by means of a metric 
tensor we also obtain a notion of density of channels. 
Such density thus provides a measure of how many dif- 
ferent channels can be distinguished in a neighborhood 
of a point S(X) with a given amount of resources 4>*. 
Such distinguishability-derived density is captured by the 
Jeffrey's prior [22j in the case of probability distribu- 
tions and the Bures prior [l4|, tDJ f° r the set of quantum 
states. These densities play a central role in the theory 
of Bayesian estimation [23[ . 

Given two channels Si and £2 and a specific probe 
state po, the maximal statistical distance [33| attain- 
able by any quantum measurement between the states 
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resulting from the action of the channels, p = S\Pq and 
a = S2P0 is given by the Bures distance [lj, [l4| 

d%(p,<r) = 2(l- y/F{p,a)), (2) 

where 

F(p, a) =tiyj y/pay/p (3) 

is the quantum Fidelity. For infinitesimally close chan- 
nels, Si — S(X) and S2 = S(X + dX) the Bures distance 
can be expressed as, 

d 2 B (p,a) = ^J(X\p ) fiU dX^dX v (4) 

where J'(X\po) is the quantum Fisher information ma- 
trix (QFI) of the model \p{X) = S(X)p } ^MM,^- 
The QFI matrix can be computed from the symmetric 
logarithmic derivatives (SLD) A^, which are the Hermi- 
tian operators that satisfy the equation 

d„p{X) = A M (X) o p(X), (5) 

where we have introduced the symmetric product for op- 
erators, A o B = \{AB + BA). The QFI is then 

J^{X) =tr[pA M oA„], (6) 

where we have dropped the explicit X dependency in p 
and A^. The SLD's play an important role in the theory 
of quantum inference |19| . 

The QFI has been reviewed a number of times in the 
literature, from which we emphasize [n], [lU, [34| [3f| . It 
enjoys several useful properties 

• The SLD has zero expectation, 

tr[pA M ] = tr[p o A„] = d^txp = 0. (7) 

• The QFI is real, symmetric and positive semi- 
definite, 

e"j flv (x)e v = tr[ /9 (rA^) 2 ] >o (8) 

which follows from p > and (# M A M ) 2 > 0. 

• The SLD and the QFI are covariant quantities. 
Given a parameter 0(X) such that dgp(X) = 
9 fi d fl p(X) we can define the SLD associated to 8, 
Ag = 6^ Afj,. We can equally define the QFI associ- 
ated to 6 as J e = 6»J^ u 6 y . 

• The QFI is monotonic under completely-positive 
trace-preserving (CPTP) maps. Writing J{p)^ y = 
tr[p A£ o A£] we have 

J{p) > J(£p), (9) 

where £ is any CPTP map. In particular, if £ is a 
unitary map, it holds that J{p) = J{£p). 

Apart from these properties, the QFI plays a central 
role in the theory of quantum statistical inference by 
placing a lower bound on the attainable variance of any 
unbiased estimator of the parameters X [l5|, [l|| [ltl Q, 



III. A CHANNEL METRIC 

In this section we define our proposed channel metric 
and study its properties. Analogously to Eq. (U]), we ex- 
pect that the distance between neighboring channels can 
be expressed as in Eq. (TTJ). This is a requirement that is 
not met by all distance measures, and it is yet unclear 
that a straightforward extrapolation of the statistical dis- 
tance applied to channels will fulfill such expectation. 
Thus, instead of focusing on maximizing the statistical 
distance, we will search for a metric tensor that provides 
an upper bound to it. Namely, 

1) -3(A) should provide an upper bound to the Bu- 
res distance between states attainable within some 
resource budget </>*. 

Following these considerations, and sticking to our defi- 
nition of p and (j, we impose the following condition on 

Z(X) 

d 2 B {p,a)<- A Z{X) tlv dX^dX v (10) 

for any possible dX and any chosen probe state po ful- 
filling 4>(po) < (j>*. Imposing this for any dX implies that 

J{X\ Po ) <3(X) Vp s.t. 0(p o ) < 4>*. (11) 

Notice, however, that this condition is not sufficient to 
uniquely specify Z{X). In general, one would expect that 
the distance established between two points is not unnec- 
essarily large. However, a tight bound for ds(X, X + dX) 
for any dX may not be expressible the form of Eq. (fTJ) . 
This is because minimizing the distance between a pair of 
points S{X) and S(X +dX) may not automatically mini- 
mize the distance between S(X) and another S(X+dX'). 
This can be formalized as follows. Let C be the set of 
achievable QFI matrices for a given channel and resource 
budget <f>* . There may exist several matrices which are 
tight upper bounds to C, namely, that they are upper 
bounds, and no smaller matrix exists which is also an 
upper bound. This is a direct consequence of the partial- 
ordered nature of matrices. There is, however, a natural 
way to precisely specify Z(X) while respecting all the 
above criteria. One expects that 

2) The volume element specified by Z(X) should be 
minimal. 

This translates into imposing that \Jdct$(X) is minimal. 

Thus, our metric corresponds to an upper bound on the 
attainable QFI for any set of Gaussian quantum channels 
£, tested with Gaussian probe states fulfilling some regu- 
larization condition 4>(po) < (j)* , and providing the small- 
est possible volume element in the manifold of channels. 
Restricting to the Gaussian domain yields a manageable 
parametrization of the probe states and allows to obtain 
explicit formulae for the QFI. In addition, the geometry 
induced by our metric will have an immediate practical 
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interpretation. It is a bound to the cfistinguishability 
attainable by Gaussian states. Moreover, recent anal- 
ysis suggests that entangled Gaussian states arc as 
sensitive as single-mode de-Gaussified states [4l[ in quan- 
tum statistical inference. This claim is yet to be proven. 

As a first approach to the problem of defining a channel 
metric it is relevant to discuss the Jamiolkowski distance 
introduced in l8| ■ It is well known that the Jamiolkowski 
isomorphism [241 ] can be used to encode a (i-dimensional 
quantum channel S into a e? 2 -dimcnsional quantum state 
ag = {S (£> 1)$, where $ is a maximally entangled state. 
One may consider that distinguishability between as 1 
and <7s 2 be related to the distinguishability between Si 
and 52. The state however, is not likely to comply 
with any physically motivated resource restriction 0. In- 
stead, one must consider the attainable distinguishability 
within the limits imposed by <j>. 

We proceed, in the next subsections, to formally de- 
fine the metric tensor Z(X) based on the previously an- 
nounced criteria 1) & 2), and prove the covariance, sta- 
bility and computability properties. 

A. Defining 

Given a parameterized family of 7i-mode channels £ = 
{5(A)}, let us define the quantum model Ai[£\p] — 
{p(X) = S(X)po}, i.e. the parameterized set of all pos- 
sible quantum states resulting from an initial probe po, 
under the action of the set of channels <£.. The associated 
QFI is 

J^(S(X)\p Q ) = tr[p(X) A p o A,], (12) 

where A M are the associated SLD's. This corresponds to 
the situation depicted in Fig. [TJa. 



a), 




FIG. 1: (Color online) Measurement schemes a) without an- 
cillas and b) with ancillas, combined with the most general 
possible measurement. Scheme b) includes a) and hence, must 
be equally or more efficient than a). 



Consider next the m-completed channel as S* m (X) = 
S(X) ®T m where T m represents the identity channel on 
m ancillary modes. We thus obtain the m-completed set 
€* m = {S* m (X)}. Given an (n + m)-mode probe state 
p we can define a new quantum model, A / f[£* m |/9o] = 
{p* m {X) = S* m {X)p }. The QFI for the new model is 

J^(S* m (X)\ P() ) = tv[p* m (X) A* m ° A*/"], (13) 

where p* m (X) = S* m (X)po is the probe state under the 
action of the m-completed channel, and the SLD's A* m 
obviously correspond to the model .M[£* m |po]- The ex- 
tension of £ to £* m is depicted in Fig. [Tf). 

Now we turn to the constraint <f>. Let us denote the 
n modes upon which the channels 5(A) act as ak,k G 
{l,...,n}, and collectively referred to as A while the 
ancillary modes are denoted by bk,k G {1, . . . , m}, col- 
lectively denoted B. We may find it useful to consider yet 
another set of modes, collectively denoted as C. We de- 
fine the acceptable constraints, as those that only involve 
the reduced density operator in modes A, pa = tr_^p, 
namely 

<i>{p) = 4>{ PA ). (14) 

This choice is physically motivated by the fact that rele- 
vant resources should only involve degrees of freedom ac- 
cessible to A. A natural choice for the function <fi could be 
the average photon number in A, <j>(p) — tr[p ^ fe aj.afc], 
which has been used previously as a comparison refer- 
ence for several estimation problems. Nevertheless, our 
results are general and not restricted to this particular 
choice. 

Finally let C(S(X)\tf>*) = {J(S(X)\p) \ <j){p) < 0*} be 
the set of all QFI's achievable by any Gaussian state ful- 
filling the constraint cf>(p) < <j>* , and C purc (S(X)\<j>*) = 
{^7(5(^)1^) | (j)(p) < (j)*,p 2 — p} the corresponding set 
when restricting to pure states. Let C be any set of pos- 
itive semidefinite bounded matrices, and J G C ' . Define 
the function M(J) = {j > J} to be the set of all upper 
bounds to J, and M(C) = n JeC M(J) be the set of all 
upper bounds common to all matrices in C. 

We are now in the position to introduce the metric for 
the set £. Let S be a channel in C The metric at point 
S is 

Z(S) = arginf det j. (15) 

{jGA/(C(S*™|0*))} 

That is, the matrix with smallest determinant that is 
greater or equal to all possible QFI's achievable by any 
2n-mode Gaussian probe state fulfilling the constraint 
4>{p) < </>*■ This is a positive semidefinite matrix and 
certainly qualifies as a metric. It is measurable, in the 
sense that it is a function of the parameters X which 
are themselves measurable [see Sec. IIV Aj . Notice that 
in very pathological cases, there may be more than one 
solution to the minimization problem. These are highly 
symmetric and unlikely situations which we will not deal 
with. We dedicate the following subsections to discuss 
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the most relevant properties of 3, and provide an ap- 
proximation method to compute it. 

The idea behind this definition is to provide an up- 
per bound to the achievable Fisher information for any 
parameter of interest, under the constraint 4>{p) < cj>*. 
While any matrix greater than those in C{S* n \(j)*) would 
certainly qualify as an upper bound, the partial-ordered 
nature of matrices prevents from having, in the general 
case, a well defined suprcmum. Instead, the minimiza- 
tion of the determinant in Eq. (|15[) aims at reducing the 
volume element to a minimum while still providing an 
upper bound to the attainable Fisher information. 

B. Covariance 

Notice that under a reparametrization X — > X'(X) we 
have 

A M A M , = A^,, (16) 

where A M ^, = dX 11 jdX'^ , where primed indices refer 
to the new coordinates and unprimed indices correspond 
to the old coordinates, as is customary in the notation 
of General Relativity. Eq. (|T5)) in turn implies that the 
QFI is a covariant quantity 

^^j; v = ^a;-av, (17) 

Let us denote this transformation law in shorthand no- 
tation J A T JA. All elements in C(S* n (X)\<t>*) 
transform in the same way under a given reparametriza- 
tion. Accordingly, the set M(C(S* n (X)\(j)*)) also 
transforms covariantly. Namely, for any matrix j £ 
M(C{S* n {X)\<t>*)) we have 

A T jA e M(C(S* n (X')\(j)*)), (18) 

and viceversa. 

Finally the determinants of all elements in 
M{C(S(X')\(f>*)) are just those of M(C(S{X)\$*)) 
multiplied by a factor detA T A. Therefore, assuming 
that under parametrization X the minimum was achieved 
for 3, the new minimum under parametrization X' will 
be achieved by 3' = A T 3A [3J,,„, = V A V A ^']' as 
required. 

C. Stability 

The stability requirement has been discussed previ- 
ously in the literature d, Q. We dedicate this section 
to prove that 3(5) is invariant under the addition of an- 
cillas, namely 

3(5 ®2) = 3(5), (19) 

where 2 is the identity channel on an arbitrary number n! 
of modes. We remark that while adding ancillary modes, 



we still require the constraints <f> to act on the original 
set of modes. Observe that 

3(5® 2)= arginf detj. (20) 

{ieM(c(S*(»+=»')|0*))} 

Observe that if X C Y then M(Y) C M(X). Given 
two sets of positive semidefinite matrices X and Y, we 
say X >z Y if and only if for any \/y £ Y 3x £ X such 
that x > y. This binary relation is reflexive (X >z X) 
and transitive (X y Y and Y h Z => X h Z). Also 
X t. Y => M(X) C M(Y). We will say that two sets 
arc equivalent (X ~ Y) if and only if M(X) = M(Y). 

Let n be the number of modes upon which 5 acts. The 
proof consists in two main steps. 

1. We will first show that the set of QFI matrices at- 
tainable using pure probe states with m > n an- 
cillary modes is equivalent to those attainable with 
pure probes using only n ancillary modes. 

2. We will then show that the set of QFI matrices 
attainable with (n + m)-mode states (m > n) is 
equivalent to those attainable with pure (n + m)- 
mode states. 

Proof of 1 : Consider an (n + m)-mode setup such as the 
one in Fig. \Tb ( m > n). The Schmidt decomposition- 
like theorem p8M51| states that, given an (n + m)-mode 
pure Gaussian state \4>}abc separated into parties A (n 
modes through the channel), B (n modes) and C (m — n 
modes), one can always reduce the system to n pairwise 
squeezed states between A and B, and m — n vacuum 
states in C by means of local symplectic transformations 
on A (5a) and joint symplectic transformations on B and 
C (5gc) [Fig-UJa]. Therefore, for any pure Gaussian state 
\?Jj)abc there is another 2n-modc state |V'2n)xB such that 

S BC \i'}ABC = \lp2n)AB\0)c (n - m) . (21) 

Defining £ as £p = SjgcpSg AC and noting that 

S <n+ m ) gp = £ S *{n+m) p WQ havCj foUowing from thc 

unitary invariance of the QFI, 

J(S* m \\i>)ABc) - J(S* m \£\i>) A Bc). (22) 

Moreover, using S* m {p (8 |0) {Qf>( m ' n )) = (S* n p) ® 
|0)(0|«'( m -™) together with Eq. (ETJ) and the additivity 
of the QFI [obviously J r (2||O) 0(m " ,l) ) = 0] we have 

J(S* m \£\ij) AB c) = J(S* m \\ij 2n }®\0)^ n -^) 

= J(S* n \\^ 2n )AB), (23) 

while 4>(\ip)) — 4>{\4'2n))- This shows that 
C pmc (S* m \cj)*) C C purc (S* n \(j)*), which combined with 
the trivial inverse inclusion yields 

C purc (5' tm |r) = C purc (5^|r), Vm>n (24) 
and therefore 

M(C pure (5* n |<^)) = M(C puTe (S* m \<l>*)), Vm>n ■ 

(25) 
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This makes it clear that all C pure (<S* m |<^*) with m > n 
are equivalent. This result has a similar flavor to the well 
known fact that maximally entangled states are optimal 
for estimation of <i-dimensional unitary operations and 
several other finite-dimensional channels [52[, which 
means that a d-dimensional ancilla is sufficient. In our 
case, maximal entanglement will in general be forbidden 
by the resource constraint <f>*, but nevertheless the 
ancilla need not be larger than the system itself. 

Proof of 2: Notice that for any (n + m)-mode QFI 
J{S* m \p) (to > n) there exists an (n+m+m')-mode pure 
state |V>> such that J(S*( m+m ">\\ip)) > J(S* m \p). To see 
this, one simply needs to construct a pure (n + m + m')- 
mode Gaussian state such that trc (V^l = Pi with C 
denoting the purifying extra (to') modes [Fig. O*]. Since 
S* m p = tr c S* (m+m 'ljiO('i/'| the monotonicity of the QFI 



under CPTP maps [53( guarantees that 

J{S* {m+m ")\\i>)) > J{S* m \ P ). 



(26) 



Notice that the constraint <p(p) is not affected by this 
construction, as long as it depends only on the reduced 
state pa = tr^p, since <fi(p) = (t>(\ip)). This shows that 
the optimal probe states can always be taken to be pure, 
provided that one enlarges sufficiently the set of ancil- 
lary modes. Therefore C pul0 (S< m+m '*> |0*) h C{S* m \<f)*) 
yielding M (C pul . c (S* ( - m+m "> \</>*)) C M(C{S* m \(j)*)). 



Using Eq. ([23)) we can, furthermore, say that 
M(C pnle {S* m \cf>*)) C M{C(S* m \cj)*)). On the other 
hand, it is trivial that C pmc (S* m |0*) C C(S* m \(j)*), 
yielding Af(C(5* m |0*)) C M(C pure (S^l <f>*)). Conse- 
quently 



M(C puie (S* m \^)) = M(C(S* m \<P*))- 



(27) 



This concludes the proof. 



This result is also a Gaussian version of previously 
known results for the estimation of finite-dimensional 
channels [54|. The added value of this result is that 
the restriction to pure states can be made preserving the 
Gaussian character of the probe states as well as respect- 
ing the resource budget (j>*. 

All these results imply that M(C{S^ n+2n '^ | </)*)) = 
M(C pure (<S*("+ 2 "')|0*)) = M(C purc (S*™|0*)), so that 



a) 



S* n {X) 




a) 



b) 




b) 




FIG. 2: (Color online) a) the Schmidt decomposition theo- 
rem ensures that in an [n + m)-mode state (m > n), one 
can locally bring m — n modes to the vacuum. Therefore, 
for any pair (\ip),Ag) of n + m modes there exists a pair 
(Sbc IVOj S^AbSbc) of 2n-modes which performs equally well. 
The remaining modes can be traced out of Sbc\iP) leading 
to a generalized measurement, which can be described by a 
POVM. The optimal POVM on the 2n-mode state is given 
by the corresponding SLD, Ag n . 



c) 




FIG. 3: (Color online) Measurement schemes a) with a ther- 
mal state in a 2n-mode setup, b) Equivalent setup with a 
pure state containing ancillary modes that are not measured, 
c) The optimal measurement including the ancillary modes. 
Obviously, the setups in s) and b) are equivalent, while the 
setup in c) has more freedom and includes b). Consequently 

•J a — Jb ^ Jc- 
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Eq. (HO]) reduces to 
3(S ®1) = 



arg inf 



detj=3(S). (28) 



b'eM(c purc (5*"|0*))} 

Summarizing, we have shown that any Gaussian state 
of 2n modes (or more) can be reduced to a 2n-mode 
Gaussian pure state which performs equally well or bet- 
ter than the original. The stability property thus follows. 



D. Computing 3(X) 

We now turn to the problem of computing the metric 
3{X). Take S as the channel for which the metric needs 
to be computed, under the constraint <fi < <f>*. Let C* = 
C(S* n \cj)*). The problem can be written as 



minimize det j 
subject to to > J yj eC*. 



(29) 
(30) 



This can be cast as a semi-infinite programming prob- 
lem (55| . where there are finitely many variables but 
infinitely many inequalities to satisfy. In the following 
we provide an approximation method by discretizing the 
problem into a convergent sequence of convex program- 
ming problems. Despite the apparently untractable na- 
ture of the problem, consistent approximations can be 
computed from the following prescription: Generate a fi- 
nite subset of constraints by sampling a subset C n C C* 
of n random QFI matrices from C* . Let Zn be the so- 
lution to the discretized problem with matrices in C n . 
Then Zn can be computed following standard convex op- 
timization methods. In fact, the problem of computing 
Zn can be recast as the problem of fin ding the Lowner- 
Jones ellipsoid of a union of ellipsoids [56| (see Fig. Q|, 
for which efficient methods exist. 

We show in Appendix [X] that this approximation 
method converges to the true value 3- More precisely, 
we show that for any e > there exists a sufficiently 
large n such that Pr(||3„ — CJ|| > e) < fcexp(— n), where 
|| • || is the operator norm. 

This iterative method assumes that one can generate 
any QFI matrix in C* . This is indeed a nontrivial task. 
In principle, QFI matrices can always be computed nu- 
merically for finite dimensional quantum systems, but in 
the case of infinite dimensions, a general method does 
not exist. In the following section we concentrate on a 
particular class of channels and provide closed analytic 
formulae for the QFI matrix. 



IV. EXPLICIT FORMULAE FOR DISSIPATIVE 
CHANNELS 

We dedicate this section to the issue of measurabil- 
ity and computability stated in the introduction. In 
principle the SLD is defined for any possible model 




(AyF 



FIG. 4: (Color online) Illustration of the problem for the 
model X = (7, N, 0, 0) at the point (7, AT) = (0.1,1), tested 
with Gaussian states fulfilling 4>(p) = tr[a'op] = 0.2. Ellipses 
correspond to the points XJ~ X X = 1. The available energy 
is distributed among single mode squeezing (red), two-mode 
squeezing (green) and displacement (blue). Ellipsoids corre- 
sponding to the QFI's for different probe states are drawn 
with matching RGB color scheme. The thick-dashed ellipse 
corresponds to the metric 0- 



S(X)p, and thus the parameters X are always measur- 
able. Similarly, the QFI is always defined and can be ob- 
tained from the second derivatives of the quantum fidelity 
F(p(X),p(X')) w.r.t. X' at X' = X 0. Therefore, the 
matrix J{X) is always defined. However, computing the 
QFI in the general case can be extremely challenging. 
Generic formulae can be found in the literature, but ob- 
taining explicit expressions often requires knowledge of 
the eigenbasis of S{X)p. On the other hand, brute force 
evaluation of the SLD amounts to solving a Sylvester 
equation [57] . Several algorithms exist for performing 
such a task numerically. Unfortunately, this is of little 
use in the case of infinite dimensions. 

In this section we develop explicit formulae for the 
computation of the SLD and the QFI of dissipative chan- 
nels probed by arbitrary Gaussian probe states. Dissipa- 
tive Gaussian channels arc of the form 



S{X) =cxp£(X), 



(31) 



where 



fc=i 



+M£D[a k }+M k D[al 



(32) 



with L[o]p = 2opo i 



op 



po*o and D[o] = 2opo 



o 2 p — po 2 and a k , a k are creation and annihilation oper- 
ators in mode k, satisfying the commutation relations 
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,ti 



Sij. We have used the shorthand notation 



X = ©Li-^fc = (7i,JVi,»A/i,9Mi, 7 2, ...) while x k 
stands for x k = (0, . . . , 0,7*, N kt &M k , SM k , 0, . . . , 0), 
with nonzero parameter values only in mode k. We have 
chosen to work with (SftM, 3M) rather than (M, M*) to 
ensure hcrmiticity of the SLD. These are the most gen- 
eral Gaussian dissipativc channels, where the {7, . . . , 7fc} 
parameters can be regarded as coupling strengths, TV cor- 
responds to the mean photon number of the reservoir 
modes and the M's are their squeezing parameters. No- 
tice that this parameterization can be taken beyond the 
natural interpretation of the parameters, e.g., by setting 
7 = c £— s/ and N — srf /{^—sf) one can account for am- 
plification processes where 3$ and 'H> are the gain, sat- 
uration and decay rate respectively in the linear regime 
{SB = 0) Moreover, the squeezing parameters M 

can accommodate phase sensitivity in the amplification 
process. 

The dynamics of an arbitrary state p undergoing the 
most general Gaussian dissipativc evolution, in the inter- 
action picture and within the Markovian approximation, 
can be described by the master equation [59[ 



(33) 



and evolves, after a time t, from state po to p = 
exp(tQ(X)) po- For all practical purposes the coupling 
constants always appear as 7fct, thus the time t can be 
absorbed in the coupling constants 7^, allowing to take 
t = 1 without loss of generality. 

There are several ways to compute the QFI. As said 
above, one consists in computing the Hessian of the 
fidelity This, however, does not provide information 
about the optimal observables, namely, the SLD's A M 's. 
Hence, we take here the longer route, by first computing 
the SLD's. This settles the issue about measurability, as 
we explicitly provide the optimal observables to estimate 
the value of X 1 *. Then the QFI is immediately given by 
their covariance matrix. In order to preserve the clar- 
ity, we introduce notation by reviewing some well-known 
facts about Gaussian states. 

At variance with the previous sections, where we have 
made explicit the distinction between the channel modes 
A and the ancillary modes B,C, we now treat all modes 
indistinctly. Considering channels with ancillas only 
amounts to considering parameter spaces of the form 



X 



ZfcQO, 

k=l 



(34) 



where = Q)™^ (0, 0, 0, 0). We thus we have n bosonic 
modes with annihilation operators a k fulfilling the canon- 
ical commutation relations [<2j,eiJ] = 5ij, all other com- 
mutators being zero. We arrange all operators into a 
vector x = ( a ij a \i ^2) c4> • • •)• The commutation rela- 
tions are expressed as 



where Q, = Q)' k=1 w, and ui is the symplectic matrix 



1 

-1 



-uj, uj 2 = -1. 



(36) 



Equivalcntly, the canonical operators Q k and P k are ar- 
ranged in the vector R = (Qi, Pi, Q2, P2, ■ ■ ■), related to 
X by R 1 = H l jX 3 where H is the unitary transformation 



k=l 



1 1 

-i i 



(37) 



and with commutation relations = iST y = 

//',.// Mr 

A Gaussian state is defined as a state with Gaussian 
characteristic function [60j . which is fully described by 
the first and second statistical moments, 

(R l ) = tv[pR% (38a) 
S lJ = tY[p{R l -{R l ))o{W (38b) 

Williamson's theorem [6l| ensures that a Gaussian state 
can always be expressed as a thermal state under the ac- 
tion of a symplectic transformation S and a displacement 
operator D, 




with 



Pv k = z k 1 CX P \-Pk a{a k 



(39) 



(40) 



where Z k is a normalization factor, and the (3 k are in- 
verse temperatures. The thermal state has zero first 
moments (R 1 ) 



0, and a covariance matrix 



1 n 

s th = ■= ffcla, v k = coth (/3 fc /2) 



(41) 



k=l 



In order to deal with states with nonzero first moments, 
it is convenient to define the displaced bosonic operators 
X 1 = Dx l D^ = x l — (x 1 )- The displaced canonical oper- 
ators R 1 arc defined likewise. With this, the CM reads 



= trlp&oRi}. 



(42) 



Symplectic transformations are linear in the cre- 
ation/annihilation operators and preserve the commuta- 
tion relations, 



S^x'S 



u l 3 x\ 



(43a) 
(43b) 



where it is a matrix representation of the symplectic 
transformation S. The latter equation is often written 
as 



[x\x J 



(35) 



u Q u T = fl, 



(44) 
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and wc will use the fact that 



u T n = flu- 1 . 



(45) 



In analogy with (|43ap , the canonical operators transform 
according to 



(46) 



where we have defined s = HuH^ , fullfilling the same 
relations (01]) and (gSJ). From (g2J| and (|35j). the trans- 
formation rule for the CM under a symplectic transfor- 
mation p U pU^ takes the explicit form 



Due to the antisymmetry of it is important to be con- 
sistent in the way that indices are lowered and risen. In- 
dices are lowered by contracting with the first index of 
Clij and risen by contracting with the second index of 
fitf, so that = ± l k W k = Y, a Vli k VLi k = E 4 '<5/. 

Finally, given a probe state po, the action of the chan- 
nel is 

p(X)=S(X)p , S(A) = exp£(A), (56) 

and the derivatives of p w.r.t. the channel parameters X M 
can be neatly expressed as 



E 



u 
u 



sis-ivE^" 



(47a) 
(47b) 



and from Eqs. (|39p and (|47a[) we have £ = s £ t h s T . 

It is worth stressing at this point the doubly- 
contravariant character of the covariance matrix £. Al- 
though the matrix equation (|47ap suggests otherwise, it 
is clear from Eq. (|47b|) that both indices transform in 
the same way. In order to construct functions of the co- 
variance matrices with proper transformation rules, we 
will find it suitable to define the covariant-contravariant 
version, by lowering one index with the symplectic ma- 
trix [48J. Let us introduce the symplectic matrix with 



lower indices fL 



rr j , so that n^n 



the 1-1 covariance tensor as 



E = Ef2, or 



E tfe 



-5\. Define 



ki- 



rn 



With this definition we see that, at variance with 
Eqs. (|47a[) and (|47b|) . under a symplectic transformation, 



s E s 



(49) 



Products of 1-1 tensors enjoy proper transformation 
rules, 

[AB}\ = A\B\, (50) 

s(AB) s - 1 = (sAs'^isBs- 1 ). (51) 



Tensor products A (g> B are 2-2 tensors with indices 



[A®B] 1 



kl 



(52) 



Moreover powers of a k-k tensor transform in the same 
way as T, 



[rpnyi . . .ik 

(s® fc T(s- 1 )® fe )™ 



j>H...l k rpr 1 ...r h 
r'i-r' k r'{...r[ 



■ T , 



The identity k-k tensor is l m , with [I] 1 j = S' 1 
the inverse of a tensor T is defined as 



r ("-i) 
k (53) 

(54) 



Finally 



T y 



(55) 



a„p(X) = v, P (x), 



(57) 



where 2? M are superoperators whose expressions are given 
in the Appendix [B] 



A. Symmetric Logarithmic Derivatives 

From now on we consistently drop the dependency on 
X and assume throughout that we are considering a par- 
ticular point S(X) £ £. The common structure of the D 
superoperators, 



Vpp = owj(xVx j - ix 3 x % ) ° p)- 



(58) 



comes in very handy for computing the SLD's in a gen- 
eral manner. From Eqs. (O and ([57]) . these satisfy the 
equation, 

V^p = A M o p. (59) 
Using x'X* — X 1 ° X* ~ we have 

V^p = a^ij X l PX J - [x l o x 3 - o p . (60) 

Combining this with Eq. (|59[) we get 



X 1 PX 3 = U M + a Mi « (x l o X J - -fF ) ° p. (61) 



This is a particular form of the Sylvester equation [5] 
Y = Z o p, which, for p > has the formal solution 



Z = 2 



'-""Ye -1 " dv. 



(62) 



Therefore 



A, 



(63) 



2 / e- v » X l PX J e- vp dv - (x l ° ^ - |n«) 
o 
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We show in Appendix [Cl that 



e- up x i pX j e~ up du 



Al X 



(/(E)® a + a<g> /(S)- 1 )- 

+ [(a + /(s))- 1 ]V J R 4 "( J R J "> 
+ [(a + /(s)- 1 )- 1 ] / ) „( J R 1 ')^'' 



R k oR l + -n kl 

ki \ 2 



+ (64) 



where E = E17 is the 1-1 covariance tensor of p and 

Finally, defining 

5^=^^®^]% (66) 
allows to express the SLD as 

+ [L^Y 3 kl R k oR l ), (67) 

where 



We conclude this subsection by commenting on possi- 
ble difficulties when computing SLD's for singular states. 
This situation can arise when the CM of p in Williamson 
form contains a vacuum mode. This is the case, for in- 
stance, when probing a zero-temperature channel (N = 
0) with a two-mode squeezed vacuum. In general the 
SLD is not defined on the kernel of p, i.e., PAP is unde- 
termined by Eq. ([5]), where P is the projector on kerp. 
Let V[A] = A — PAP. By measuring A v on the state 
p(X + dX) = p(X) + p(X) o A^dX* it is easy to see 
that tr[p(X + dX)A v ] = tr[p(X + dX)V[A v }\. This is 
the freedom available in defining SLD's for singular p's. 
However, Eqs. (frJT)) and (j6"5| were derived from Eq. ([55]) 
by means of expression (|6"2"|) which assumes that the den- 
sity operator is nonsingular. Thus, these expressions can 
yield to divergencies for PA^P. These are not observ- 
able divergencies, and can be regularized by introducing 
a small temperature e in the probe state po, projecting A 
with V , and finally taking the limit e to zero. There may 
be, on the other hand, observable divergencies. This sit- 
uation often arises when one explores the vicinity of the 
boundary of the manifold. The classical Fisher informa- 
tion is well known to diverge in many statistical models, 
the most prominent case being the binomial distribution 
when p — >• or 1, giving rise to Poissonian statistics. 
Analogously, the Burcs distance for mixed qubit states is 
well known to diverge in the limit of pure states, giving 
rise to interesting and counterintuitive effects in Bayesian 
qubit estimation (63^ . A similar effect in quantum statis- 
tics has been reported in a particular Gaussian channel 
[X = (7, 0, 0, 0)], when the parameter T ->■ 



=i([f(t) <g> a + a <g> /(— E)] _1 + -a <g> a), 

(68a) 

[LW]>[2(l + /(E))- 1 -l] i k (ff) 

+ [2(1 + /(E)- 1 )- 1 -1]^), (68b) 

l( 2 ) =2([/(e) <g> a + a <g> f(-t)}- 1 - ^1 <g> 1). 

(68c) 

Here we have used the fact that f(x)^ 1 = f(—x). 

Notice that this expression overcomes the main dif- 
ficulty in evaluating the QFI for Gaussian channels of 
continuous variable systems, namely, expressing the SLD 
in a manageable form. In our case, the problem is re- 
duced from working in an infinite-dimensional Hilbert 
space, to dealing with finite dimensional vector spaces, 
where the tensors L^ ', ZA 1 ) and arc defined. More 
importantly, these tensors depend exclusively on the co- 
variance matrix and the first moments of the quantum 
state output from the channel, which have a simple re- 
lation to the channel parameters [62j . This will not only 
provide a means for evaluating the channel metric 2, but 
it will also allow to evaluate the performance of several 
channel measurement schemes 1401. 



B. The Quantum Fisher Information matrix 

The expression of the SLD, Eq. ([BT]) . allows for the 
computation of the QFI, Eq. ([6]). The full derivation of 
the most general QFI is given in Appendix [El Eq. (|E13[) . 
We reproduce here the resulting expression: 

J^{X)=a^ f a v , kn , ([£ (1 >]*'{ [L^] k '{^ k (69) 

+ [l (2) ]%„ [L {2) ]^[Df-f {a lk w l + n«n*) Y 

It should be noted that the tensors appearing in Eq. (|6"9"|) 
depend on E, the covariance tensor of the output state. 
In order to finally evaluate the QFI one needs to express 
E in terms of the initial quantum state Eo and the chan- 
nel parameters X. The corresponding relation is well 
known in the literature [62l [65[ but we repeat it here for 
completeness. Arranging the channel parameters X in 
the asymptotic covariance matrix E c h, 

v - n\( \+ N k + ®M k QM k \ , 
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and defining the coupling matrix 

r = e -^l 2l (71) 
k 

we have 

E = r(So-s C h)r + s ch , (72) 

where So is the CM of the initial state, and £ is the CM 
of the state going out of the channel. 

Eqs. ([11]), (HOI), CLE]) and ([72D provide a systematic way 
to compute the QFI for any channel S(X) and any input 
state po- Remarkably, this expression is exact and an- 
alytical. However, despite the significant simplification 
of the problem - namely, from an infinite dimensional 
operator equation (|59p to a finite-dimensional matrix ex- 
pression, Eq. (|69[) - explicit analytical expressions are too 
complex to be of any use, except for very simple channel 
models. 



V. CONCLUSION AND OUTLOOK 

Summarizing, for any resource constraint (j){PA) < </>* 
we have defined a metric tensor on the manifold of Gaus- 
sian channels. We have proven that the metric is stable 
under the addition of ancillas and provided a method for 
computing numerical approximations, by using convex 
optimization methods. The resulting distance is an upper 
bound to the attainable Bures distance between states 
resulting from the action of the channel onto any given 
initial Gaussian state fulfilling the constraints, and it al- 
lows to establish a systematic way to measure distances 
between channels in order to quantify imperfections in 
quantum information implementations. Moreover, the 
metric tensor minimizes the volume element assigned at 
any point of the set of Gaussian channels. This density 
can be used as an analog of the Jeffrey's (Bures) prior 
distribution for points in the simplex (density matrices) 
in Bayesian estimation methods. Our results are a step 
towards the identification of a useful notion of distance 
among bosonic channels, which are the basis mathemat- 
ical framework of quantum communication devices with 
continuous variables. 

Additionally, our approach has provided several results 
which are of relevance to the field of quantum estimation 
theory. We have derived closed formulae for computing 
the QFI of any estimation method that uses Gaussian 
states and dissipative Gaussian channels. Moreover, for 
each one of the channel parameters we have obtained 
closed expressions for the corresponding optimal observ- 
ables associated to each one of the parameters (symmet- 
ric logarithmic derivatives). Moreover, we have proved 
that by enlarging sufficiently the number of ancillary 
modes, the optimal probe states for estimating Gaussian 
channels are always pure. 

These findings should be of immediate use in several 
practical situations. The most prominent of them is 



entanglement-enhanced metrology of Gaussian channels. 
A first relevant application of our results is the evalua- 
tion of Gaussian protocols for parameter estimation in 
dissipative Gaussian channels, such as decay rate, tem- 
perature or degree of squeezing in engineered baths (40j . 
The stability proof in Sec. MI CI establishes the maxi- 
mum number of ancillary modes to be considered for any 
of these situations. Indeed, such proof is very general 
and can be extended to several metrology problems, lim- 
iting the number of ancillary modes required for opti- 
mally, thus restricting the number of setups that need 
to be explored when designing optimal metrology pro- 
tocols. In the spirit of the Choi-Jamiolkowski theorem, 
our arguments show that estimation of single mode chan- 
nels cannot be enhanced by multipartite entanglement as 
compared to bipartite entanglement, at least within the 
Gaussian framework. 

Moreover, the approach taken here may be extended 
to derive SLD's for other related problems such as phase 
estimation under decoherence, phase diffusion, etc, which 
can be expressed quite naturally in phase space, and that 
would otherwise be difficult to address due to the infinite 
dimensional character of continuous variable systems. 

Our work leaves open several questions and possible 
extensions. Obviously, obtaining explicit analytic ex- 
pressions for the channel metric would be the ultimate 
achievement. However, this seems to be out of reach un- 
less some significant advances are made. In particular 
we point out three relevant missing points; a) A clear 
criterion for the determination of a set of optimal probe 
states for any given choice of </>. b) A general form of 
the QFI for all Gaussian channels, not only the dissipa- 
tive ones, c) An analytic expression for the metric in 
the high resource limit (properly regularized). Concen- 
trating our attention on the regularization schemes, we 
can point out some interesting questions: What classes 
of constraints <f> will always be saturated by the opti- 
mizing states? Considering a related problem, prelimi- 
nary numerical results (40| show that in the single-mode 
bosonic lossy channel, at zero-temperature, with con- 
straints of the form (j){p) = tr[a^ap], an entangled two- 
mode squeezed state of the probe and ancilla modes with 
mean photon number (n) in the channel mode will per- 
form equally well any other state with n average photons. 
If the result of this preliminary analysis will be confirmed, 
it is unlikely to be a coincidence. Finding the reasons 
behind this surprising match is an interesting problem 
worth deeper investigation. 

Naturally, these are not the only open questions, and in 
particular, extending our work to non-Gaussian channels 
and/or non-Gaussian probes would be of utmost impor- 
tance. There is, in fact, nothing peculiar in our defining 
scheme that relies crucially on Gaussianity of the chan- 
nels and/or the probe states. The metric could thus be 
defined on more general sets of channels and probe states. 
It is however, the restriction to Gaussian channels and 
states that allows us to provide computational formu- 
lae and prove that our metric satisfies the basic require- 
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ments listed in the introduction, especially the stability 
property. This does not mean that non-Gaussian exten- 
sions will not fulfill the requirements. Indeed, it would 
be an interesting result if a counterexample to stability 
could be found in the large set of non-Gaussian states 
and channels. It would mean that some estimation pro- 
tocols can be improved by considering multipartite se- 
tups, when only one part (mode) undergoes the channel 
action. This would be an effect exclusively of continuous 
variable systems which cannot occur in finite dimensions. 

Finally, a question remains unanswered regarding 
the way in which channels are combined. Throughout 
the text we have assumed that the channel is tested a 
large number N of times, independently, namely, a given 
probe state is prepared, sent through the channel and 
measured, and the process is repeated N times. The 
number of trials is immediately regularized and thus does 
not enter in the discussion. However, in the quantum 
scenario, regularizing ./V leaves no room for chaining 
channels together, so that, for instance, a state may be 
sent twice through the channel before being measured, 
or entangled states may be used in the simultaneous 
testing of a single and a doubly chained channel. It is 
evident that a large number of degrees of freedom are 
not exploited in our approach. However, our framework 
yields the most natural bounds, which can be obtained 
in the iid ("independent and identically distributed") 
case, when no correlation is admitted between different 
samples. A very relevant and interesting open question 
is whether such sophisticated schemes can enhance the 
precision of quantum estimation and quantum metrology 
setups. 
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Appendix A: Convergence of the approximation 
algorithm 

We show here that the approximation method given in 
Sec. IIII Di does converge to the true value of the metric Z- 
Throughout this appendix we make extensive use of the 
Hausdorff distance dn (•, •) [66j|, which defines a distance 
among subsets of a metric space W. In our case W is 
the set of positive semidefinite matrices, equipped with 



the operator norm distance d(x,y) = \\x — y\\, where \\x\\ 
is the operator norm ||x|| = sup|^ =1 \xv\. Define the e- 
ball centered at a point x as the set of points within a 
distance e of x, B e (x) = {x' | d(x,x') < e}. Then, the 
e-ncighborhood of a subset X is the union of all e-balls of 
X, B e (X) = U X £xB e (x). Then, the Hausdorff distance 
is defined as 

d H {X, Y) = inf{<5 | X c B S {Y) and Y C Bg(X)}. (Al) 

Let us first show some preliminary facts about the 
randomly generated sets of constraints, C n . 

Theorem: Assuming that the set C* = C(S\<f>*) of 
QFI matrices associated to a given channel S [within 
the corresponding probe state restrictions 4>{p) < </>*] 
is bounded, we have the following property: For any 
e > the probability Pr(dy(C„, C*) > e) decreases 
exponentially. 

Proof: Since the set C* is bounded it can be covered by 
a finite number k of e/2-balls centered at points {yi € 
C*}, i = l...,k, thus C* C B t/2 {{ Vl }). If 

B e/2 ({ Vl }) c B e (C n ) (A2) 

then C* C B e/2 {{yi}) C B e {C n ). Since also C n c C* C 
B e (C*) we have that Eq. (X2j) implies d H {C n ,C*) < e. 

We now compute an upper bound to the probability 
that Eq. (|A2[) is true. To each random selection x € 
C n we assign the point y e {yi} which is closest to x. 
Since d(x,y) < e/2 then B e / 2 (y) Q B e (x). Let m be 
the number of points of C„ assigned to ball B e / 2 (yi). A 
sufficient condition for Eq. (|A2p to hold is ?i; ^ Vi 
[all points in {y} have been assigned at least one point 
x e C„]. Thus 

n^OVi => d H {C n ,C*)<e. (A3) 

Consequently, 

Pv(d H (C n ,C*) <e) > Pr(n^OVi). (A4) 

Let pi be the probability that a random sample x is as- 
signed to point yi. The probability that after n samplings 
a number m, . . . , nu of points has been assigned to each 
ball is given by the multinomial distribution 

n' 

Pr(ni "--' nfc) = ^T^! Pl - p * 

We now bound the complementary probabilities 
Vr{d H (C ni C*) > e) < Pr(3i s.t. n t = 0) 

= Pr(0, n 2 , rife) + perm. 

+ Pr(0, 0, nk) + perm. 

n 3 ,...,n k ^0 

+ ... 

+ Pr(0,0,...,0,n)+perm., (A5) 
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where perm, represents all permutations among ar- 
guments of the multinomial distribution and it is im- 
plicit that summation is over all values such that 
Sfe n fc = n - We can now upper bound Pr(3i s.t. rii = 0) 
by completing the sums to include indices equal to zero, 
and by using the multinomial theorem. 

Pr(3i s.t. m = 0) 

< (Pa + P3 + ■ ■ ■ PkT + perm. 
+ (p3 + . . . + Pk) n + perm. 
+ ... 

+ Pk + perm. 
= (1 — p\) n + perm. 

+ (1 — pi — p 2 ) n + perm. 

+ ... (A6) 

+ (1-pi -P2- ■■■ -Pk-i) n +perm. , 

To conclude we can further upper bound this quantity 
by replacing each term (•)" by the maximum value (1 — 
mim^i)" and counting the number of terms YljZi (;) = 
2 k - 2, 

Pr(3i s.t. m = 0) < (2 k - 2)(1 - minp,)™ (A7) 

i 

Notice that this is exponentially decreasing with 
n. Thus, we have that for any e > the probability 
Pr(dff(C ra , C*) > e) decreases exponentially with n. 
This means that by picking enough many samples we 
have djj(C n , C*) < e with arbitrarily high probability. ■ 

Let x € X be a positive semidefinite matrix within a 
set X. Define M(x) as the set of all upper bounds to x, 
and M(X) as the set of upper bounds to all elements in 
X, i.e., M(X) = {m\m > x G X}. We call M(X) 
the set of feasible points to the optimization problem. 

Lemma: Given two matrices x and y such that 
d(x, y) < 5 it holds that x — 61 <y < x + 51. 

Proof: Define A = y - x. Then ||A|| = d(x,y) < 5. 
Moreover -||A||1 < A < ||A||1 so that -51 < A < 61. 
Thus, 

x-61 <y <x + 61. (A8) 



Theorem: The map M : X \-t M(X) is con- 
tinuous in the Hausdorff distance. Namely, for 
any e > there exists a 6 > such that 
d H (X,Y)<6 d H (M(X),M{X)) < e. 

Proof: For any e > take S < e. We prove constructively 
by showing that if dg(X, Y) < 6 it follows that M{X) C 
B t (M (Y)) and M (Y) C B £ (M (X)) which in turn implies 
that d H (M(X),M(Y)) < e. 

We have to show that a G B e (M(Y)) Va € M(X). 
Take any element a G M(X) and construct 6 = a + 51. 



It must follow that (j) 6 G M(Y) and (ii) d(a, b) < e. 
(i) Since dji{X, Y) < 6 then for all y e 7 we have 
an x G X such that d(x, y) < 6 and by the lemma, 
x + 61 > y. Thus, b = a + 51 > x + 51 > y for all 
y G Y. Thus 6 G M(Y). (ii) d(a, b) = 6 < e. This proves 
that M(X) C B e (M(Y)). The exact same reasoning 
with X and Y interchanged shows the converse, thus 
d H (M(X),M(Y)) <e. * 

The last two theorems combined show that given 
any e > 0, one can always generate a set of feasible 
points M(C n ) which has Hausdorff distance from M(C*) 
smaller than e with arbitrarily high probability, i.e., 
Pi(d H (M(C n ),M(C*)) < e) > l-fcexp(-n). Define 
the function 

f(X) = arginf det x. (A9) 

It only remains to prove that 3« = f(M(C n )) approaches 
Z as M(G n ) approaches M(C*). Let us remark that, 
by construction M(C*) C M(C n ). As said in the 
text, we assume that the solution to the real problem 
3 = f(M(C*)) is nondegenerate. This means that 
det (J + A) > det J whenever 5 + A G M(C*) and 
l|A|| # 0. ' 

Theorem: Assuming that the miminum of 
{dcty | y G Y} is nondegenerate, the map {f(X n )}, 
converges to f(Y) in the Hausdorff distance whenever 
{X n \X„ D Y Vn} converges to Y. Namely, for any e > 
there is a 5 > such that if Y C X and dff(X, Y) < 5 
then ||/(X)-/(Y)|| <e. 

Proof: Let x* — f(X) and y* — f(Y). Consider an open 
ball of radius d < e/2 around y*, and let T be the smallest 
gap between determinants of points in Y outside Bd(y*) 
and det y*, i.e., 

T= inf doty- det w*. (A10) 

By the nondegeneracy assumption T > 0. Consider the 
function 

h(5) = n{\\y*\\+D + 5) n - 1 5 (All) 

where D is the maximum distance between points in the 
boundary of Y. Clearly h is monotonically increasing 
and h(0) = 0. Let 

5 < min(/i" 1 (r),e/2) , (A12) 

By assumption X C Bg(Y) and Y C X. Let y G Y be 
the closest point to x*. It is easy to see that y must lie 
in the boundary of Y, since if it was not, there would be 
a ball B T (y) C Y and y' = y- t\\x* - y\\l G B T (y) C Y 
would have — x*j = (1 — t)||?/ — x*|| , which contradicts 
the assumption that y is the closest point to x* in Y. 
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We know that ||x* - y| < S. Let y = y* + A. By the 
triangle inequality 

||x*-y*|| <||A||+<5. (A13) 

Also, since Y C X we have detx* < dety*, thus 

dety - dety* < dct y - detx*. (A14) 

Moreover (see (57j ) 

dety -detx* < »max(||»|| > ||x*||) n_1 ||x*-y||.(A15) 

We can further simplify this by noticing that ||x*| = 
lb + (x* - y)|| < \\y\\ + \\x* - y\\ < \\y\\ + S, thus 
max(||i4||s*||) < \\y\\ + S < ||y*|| + || A|| + S, so we have 

det y - det x* < n(\\y*\\ + \\A\\ + S^S. (A16) 

Notice that || A|| = \\y — y*\\ is the distance between two 
points in the boundary of V, thus || A|| < D, so that 



dety — detx* < h(6). 



(A17) 



On the other hand, let us show by contradiction that 
||A|| < d. Suppose y £ B d (y*). By Eq. (fATTj) . the mono- 
tonicity of h and the definition of 5 we have 



dety - detx* < h(S) < T. 



(A18) 



which means that detx* > dety — T. But y £ Bd(y*) 
means dety — dety* > T, which yields detx* > dety*. 
This contradiction shows that ||A|| < d, which in turn 
shows that ||A|| < e/2. 

Finally, recovering Eq. (|A13[) we obtain 



(A19) 



Taking = f(M(C n )) and 3 = /(M(C*)) and ap- 
plying this theorem to f (M(C n )) and /(M(C*)), where 
C n C C* and M(C*) C M(C) with 



Pr(d H {M(C n ),M(C*)) >S) 
will ensure that 







Pr(||3n-3|| >e)™0 



(A20) 
(A21) 



Since the probability to generate a C n within any given 
Hausdorff distance to C* approaches 1 exponentially, we 
see that any desired precision can be attained with arbi- 
trarily high certainty by sampling enough matrices from 
C*. 



Appendix B: Derivaties w.r.t. the channel 
parameters 

In order to compute the derivatives of p(X) = S(X)po, 
note that 



d M p(x) = dpS(x) po. 



(Bl) 



Thus we only need to compute the derivative of the su- 
peroperator S(X). Using the relation (6?J 
fl 



d 



cxpG(X) 



,uQ(X 



d li g{X)e^-^ G ^ x Uu 1 (B2) 



dX» Jo 
the channel derivatives can be written as 

d^S(X)= [ e^^d^Gix^e-^^duSiX), (B3) 
Jo 

where k stands for the mode to which the parameter 
X^ corresponds. This expression allows for the handy 
relation 



cy>(A) = v fl s(X), 



(B4) 



with 



V^= / e uS ^d„g(x K )e- ue ^du. (B5) 
Jo 

Using the Hadamard lemma, 

exp(uA)B exp(— uA) = } — -[A, B] m , 



m=0 



(B6) 



with [A,B] = B and [A,B] m = [A, [A,B] m -{\, and the 
commutation relations 



[L[a[],L[a k .\ 
[D[al),D[a k ,] 
[D[a k ],L[a k ,] 
[L[a\],D[a kl ] 
[D[a{],L[a k ,] 
[L[a{],D[al,] 



= 25 k k>{L[a k }+L[al}), (B7a) 



0, (B7b) 
28 kk ,D[a k ], (B7c) 
28 kk ,D[a k ], (B7d) 
25 w D[all (B7c) 
28 kk ,D[a\], (B7f) 
one can show by induction that 

[Q(X),L[a k ]] m = lk n L[a k ]+2 1 ™- 1 Q{x k ), {B% a ) 
[g(X),L[a{]] m = ^L[a{]-2^- 1 g{x k ),(B8b) 



[g(X),D[a k ]] m = y?D[a k ], 

[g(X),D[al]] m = ^D[al\, 

and the T> superoperators become 
1 



V 



Ik 



1>1StM k 



Ik 

(e 7fc 



-g(x k ), 



i 



(e 7fc - 1) 
(e 7& - 1) 



L[a k ] 


+ 






2 




D[a k ] 


+ 


D[a\] 




2 




D[a k ] 




D[a{] 



2 i 



(B8c) 
(B8d) 

(B9a) 
(B9b) 
(B9c) 
(B9d) 



Finally, notice that the T> superoperators allow for eval- 
uating derivatives of the state p in a simple form, 

d^p = d ti S(X) Po = V^S(X) Po = V^p, (BIO) 

where all have the structure 



T>»P = a^ijWpX 3 ~ (rV) ° P)- 



(Bll) 
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Appendix C: Computing the SLD 

We proceed here to show how Eqs. (|S7|) and 

(jnHJ) ar e derived. We remark that the relations obtained 
apply to all Gaussian states in any number of modes. 
We use the notations introduced in Sec. ITVl 



Lemma: [p th , x % ] 
fe diag(/3 fc ,-/3 fe ). 
Proof: We begin by showing 



[exp M — IJ' xVth where M 



(CI) 



Define A = - fc p k a\a k and B m = [A, B m _i]; B = X*- 
One can show by induction that B m = \M m \^ . Using 
the Hadamard lemma, 



where Co = x l an d C m = [pth,C m _i]. One can show 
by induction that C r „ = [{E - ^) m })x j Pll- Now observe 
that 



e~ vp x l e vp 



DSe~ vpth S'<x l Se vpth S'<D^ 
s\, DSe- vpth X l 'e vpth S^D^ 

(_ v \m 



z — ^ m ' 



(C6) 



pthxVu, 1 = Ei B ™ = [ ex P M ]V' ( C2 ) 



thus, finally 



Lemma: Given p > with CM E and = Dx l D\ the 
following holds 



[Pth, = (PthXVth 1 ~ X 1 ) Pth 
= [cxpM-l]^Vthl 



(C3) 



Lemma: Let u be the symplectic representation of S as 
defined in Eq. (|43a[) and s as defined in Eq. (|4l)]) . where 
S 1 is defined from p in Eq. ([31?]) . Let E be the covariance 
matrix of p (Eq. ([32])) and F = exp M. Then 



where 



uBii" 1 = H*f(£)H 



x + i/2' 
Proof: Observe that 

E 2 - 1/4- iE 



E 2 + 1/4 
E? h -l/4-iE 



th -l 
S 



e 



A- 



cosh fik i sinh 
— isinh/3fc cosh/3^ 



HuEu^H^ 



(C4) 



-^ e -^=[(l + F)- 1 ] i ^'p- 1 



^(l + F- 1 )- 1 ] 1 ^" 1 ^ (C7) 



Proof : 



e- vp x l e- vp dv 



E^ 



1 - F 



xV 1 



The other identity is proven analogously 



Theorem: Given p > 0, the following relation holds 



(C8) 



Lemma: Given p > with CM E, define F = H\f(t)H. 
The following relation holds 



e- vp x l e v 



P = E 



(-«)' 



e -« P ^ip^J e""" d« = [(F® 1 + 1® F- 1 y 1 f kl x k X \ 



Proof: Again, using the Hadamard lemma, 



where the inverse of a tensor TV ; is T 1 such that 



E 



'Co 



(C5) 
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Proof: Making extensive use of the previous lemma, 



e-^P^px 3 e~ vp dv 



L — ' ml 



x k 



x / (-v) m p m+1 e- vp x ] e~ vp dv 



E 



(-1)' 



mini 



\{F~t) m U{F-tT]\x«p 



n]j ~k 1+m+n 



v m+n e -2vp p-ntf n du 



x / v e ">■ p -'XP 
la 

(-l) m (m + n)! 



E 



2 1+m +™m!n! 



[{F-tru{F-tr]\ 



E 



(-1) TO ( m + n 



2l+m+n y 77, 

mn 

(-l) m f m + n 



E 



[(f - iry k [(F - irw 



[(F-trMi-F- 1 )™}] 



(C9) 



Next, make the change of variables m + n = q and replace 



m,n=0 ^ 2^ g =0 Z^n=0> 



= E E ( n ) - 8 c^ -1 - 

q =0 n=0 ^ ' 

= \(F ®t + l® F- 1 y 1 V J . (CIO) 

L J kl 

This proves the theorem. ■ 

Finally we are in position to obtain Eq. (f6"4"| . Observe 
that 



e vp x 1 ' PX° e vp dv 



(Cll) 



e- vp (x+(x)Tp(x+(x)ye~' vp dv. 



Collecting the results and using 

[H®H] kl m , x k 'x 1 ' =R k R l 



R k o r 1 + -o* 



[H®H] kl k >i> x k \x l ')=R k (R l ) 
[H®H] kl kll , (x k ')(x l ') = (R k )(R l ) 



(C12) 

(C13) 
(C14) 



we finally obtain 

e~ up X l PX 3 e~ up du 
[H* ® H^f vy x 

(/(E)® 1 + 1® fit)- 1 )- 1 



R k oR l + -n kl 



kl 



+ [(l + f(£))- 1 ]\ i „R*"(W) 

+ [(i + f(£)- 1 )-%,(R l ')R J ' 
+ \(R l ')(R f ?j 

which is Eq. 



(C15) 



Appendix D: Some simplifying identities 

In this section we derive some relations that will be 
useful for the simplification and evaluation of the QFI. 
Consider the tensor [/(E) ® 1 + 1 ® /(-E)] _1 . Notice 
that if [X, Y] = we can write 

lnx> + f>y)]-' = ^ XY - l/i (x + r> 



2 \XY + 1/4 217 + 1/4 

(Dl) 

Given that /(X) ® 1 = /(X <g> 1), 1 ® f(Y) = /(I ® F) 
and [E ® 1, 1 ® E] = we can write 



[/(E) ® 1 + 1® /(-E)]" 



(D2) 



1 /- i ■ 

-D _1 E®E+-1®1- -(E ® 1 - 1 ® E) 

2 V 4 2 



where we have defined 



D = E® E- -1® 1, 



(D3) 



thus 



=iL>- 1 ^E®E-^(E®l-l®E)), (D4) 

L (2) = -ZT^l® 1-«(E®1-1®E)). (D5) 

At this point, the following relations arc useful, 

[E®E] y H n kl = [EOO(EO) T ] ij 

= [EOEp, (D6a) 
[E ® 1 — 1 ® t] ij M Cl kl = [EOO - 0(E0) T P, 
= -2E lJ ' 



[E®E] l \ ; E fc ' = [EOE(EO) T ] ij ' 



[E ® 1 - 1 ® E] iJ fe; E fc ' = [EOE - E(E0) T P 



(D6b) 

[EOEOE]*- 7 , (D6c) 
:OE - E(EO) T p 
2[EOE] y . (D6d) 
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Thus finally we get 



D- 1 ( + -S 



D- 1 + -£ 



thus 



[i<°>]V w = -[^ , ]V 



(D7) 
(D8) 

(D9) 



This identity, together with the fact that tr[pR'} = 
guarantees the consistency check tr[pA M ] = 0, Eq. (JZJ. It 
also turns out to be very handy in the simplification of 
the QFI, in Appendix [El 



Appendix E: Computing the QFI 

The QFI can be computed from Eqs. ©, (p7|) and 
dHHJ) - Combining Eqs. ^ and (jUT)) we obtain all possible 
cross terms among L (0 \ and . Classifying terms 
according to the number of R operators they contain, 
we find the following. Of zero order only one term is 
obtained (|Ela[) . Terms in first order contain products of 
L(°) and ZyW. These terms vanish when the expectation 
is taken ((R) = 0). Second order terms contain cross 
products between and and "second powers" of 
L (1) , lines (|Elbj) . (|ETc|) and (|Eld|) . Third order terms 
contain expressions of the form trfpi? 1 o R )], which 
vanish identically [we argue this claim at the end of the 



appendix]. There is only one fourth order term (|Ele|l . 
which comes from double contribution of L^ 2 \ 



(Ela) 
(Elb) 
(Elc) 
(Eld) 



+ [LW}%., [L^] kl kll M[p(R 1 ' o Ri') o (R k ' o #')] 



(Ele) 



Using Eq. (|E)9|) this is reduced to 

Jnv(X) = a Ml ij& Vi ki x 

[L^] ij .,[L^] kl k ,^' k ' + [LV>]%.,1?1' [L^] kl k , v n k ' 1 ' 

+ [L^,[L^l l M[p(R'oR^o(R k 'oR^ 

(E2) 

At this point a manageable expression for T tlkl = 
tr[p(i?' o o o R 1 )] is an imperative. This ten- 
sor is symmetric under interchange of indices (12), (34) 



and simultaneous interchange (13) (24), which generate a 
subgroup (3 of the full symmetric group S±. It is conve- 
nient to establish the relation between the tensor T zjkl 
and the fully symmetric one 



ijkl 
W 



{iju) 



4! 



(E3) 



which can be easily computed with standard phase space 
methods (HI, 



T% kl =tT[pW(R'WR k R% 



(E4) 



where W(X) is the Weyl-ordered product of X. Being S 
the covariance matrix of p, T^ 1 is just the fourth mo- 
ment of the Gaussian distribution with covariance matrix 
S, namely, 



rpijkl Y^kl _j_ Y? k Yj^ -\- S^S^ 



(E5) 



In order to exploit the preexisting symmetry of T we 
split the sum over a € S4 over right cosets 54/6, of 
which we pick {e, (23), (24)} as representatives, namely 
6 U (23) e U (24)6 = S* 4 . We can rewrite T w as 



ijkl 

w 4! 



1 ^ ^ T aa(i3kl) 



(E6) 



aS{e,(23),(24)} ctGS 

Since T is invariant under the action of & we have 



ijkl 
W 



{iju) 



l©! 

4! 

ae{e,(23),(24)} 
1 (rpijkl _|_ rpikjl _|_ rpilkj^ (E7) 



Therefore 



rpijkl rpijkl _|_ ^- ppijkl pikjl^ _|_ ^ (rpijkl pilkj^ 

3 3 



rpijkl i f rpijkl , rpijkl\ 

— 1 W +3l > I [2,3] +J [2,4]J ' 



(E8) 



where 3 j is the antisymmetrized tensor in the second 
and third indices, 



T, 



ijkl 



1 



[2,3] 2 



^rpijkl _ rpikjl^ 

- (n ij n kl - a ik n jl ) - ^rt u n jk , (E9) 



and T#*l = § (T^ kl - T 1 ^') = | (T^ lk - T ? ^' fc ) 



2,4] 

T [2,3r With this 

we obtain 



kl rpijkl 



rpijkl rp 



1 



- (n ik n jl + frrp 



(E10) 



Thus finally 

rpijkl ynj y]fcfl J yjikyj^/ _|_ ^il-Y^jk 

_ (QikQjl + n ilQjk\ 

4 v y 



(Ell) 
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It will be convenient to write this result as 

where D is defined in Eq. (|D3|) . 

A similar approach shows that R 1 o (R? o R k ) = 
W(R l B?R k )+ terms containing only first powers of R. 
It is easy to see that expectation values of all these terms 
vanish identically. 

Using the expression for T l i kl we get 

+ t i(2) ] %» [ L(2) ] m'^'IT ( n<<w '' + ^ lQjk ) ) » 

(E13) 



where, by virtue of Eq. (|D9[) . the second term in Eq. (IE2|) 
has been canceled against the first contrubition coming 
from Eq. (|E12]l . 
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